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Abstract. We study the 7r ± -spectra obtained in 2, 4, 6 and 8A GeV Au-Au collisions within the thermal 
model. We find that the main features of the data can be well described after we include the pions from 
the decay of the Zi-resonance with medium mass modification. 

PACS. 25.75.-q Relativistic heavy-ion collisions - 24. 10. Pa Thermal and statistical models - 25.75.Dw 
Particle and resonance production - 24.10.NZ Hydrodynamic models 



1 Introduction 

Heavy ion collision experiments, in the energy range of 
several GeV, permit the study of nuclear matter and hadron 
properties under extreme conditions of temperature and 
density. Here, we are interested in how the pion spectra are 
influenced by the in medium mass modification. Consid- 
ering the most precise pion spectra available, , obtained in 
the energy range 2-8 A GeV [1] , we show that there is sen- 
sitivity to the medium modification of the mass splitting 
between ^33 and N. The availability of both pion charge 
polarities allows us to address the pion isospin asymmetry, 
an important ingredient in our analysis. 

A prominent feature of these 7r ± -spectra is that they 
are not described by a single thermal source. Instead, the 
initial analysis by the experimental group used a two- 
component thermal model with two different source tem- 
peratures. The dynamical origin of such a model is at 
present open to discussion. On the other hand, the de- 
cay Z\33(1232) — * N + 7r has long been recognized to be a 
significant mechanism of pion production in the reaction 
energy range considered here [2 , 3J . Moreover, the shape 
of the decay pion spectrum is different from the shape of 
a thermal spectrum [3]. 

Thus, in order to obtain the observed momentum spec- 
tra of pions, we consider the sum of the direct thermal pion 
component with the in-medium-decay of the thermal ZI33 
component. Only decays near the kinetic freeze-out con- 
tribute in the second component: the early decay product 
pions are reequilibrated and are part of the pion ther- 
mal component. The number of particles at a given value 
of momentum (the spectrum) is established by consider- 
ing statistical distributions which arise at temperature- 
dependent values of hadron masses. Below the kinetic (i.e., 
scattering) freeze-out temperature at about T = 100 ± 20 



MeV, in the free-streaming domain, the hadron mass val- 
ues are restored to the free space value. Moreover, as 
hadrons begin to free-stream out of the interaction region, 
their mass returns to the free space value, the required 
energy is derived from the modification of the volume- 
vacuum energy. Since all particles actually observed ex- 
perimentally had time to 'recover' from medium modifica- 
tions to their properties, till now the experimental study 
of this medium effect has focused on the observation of 
the possible modification of the decay width in matter, 
especially of the </>-meson [5]. 

In order to describe the pion spectra we will show that 
the required magnitude of the mass medium modification 
is in agreement with the widely held believes about the 
medium dependent properties of strongly interacting par- 
ticles in dense, hot matter [5J[7] . Considering the quantum- 
hadrodynamic model [5] , the cancellation of several contri- 
butions to the energy of a nucleon with magnitude of a few 
hundred MeV leads to the relatively small nuclear bind- 
ing energy in normal nuclear matter. One of these large 
components is the scalar potential which modifies nucleon 
mass. At normal nuclear density it has a magnitude of sev- 
eral 100 MeV, and as the nuclear matter is squeezed, this 
potential increases such that the effective nucleon mass 
melts entirely (see for example Fig 2. in Ref.(S]). Already 
at relatively moderate temperatures the thermal melting 
is more important compared to the density effect. A way 
to express the hybrid temperature -density dependence of 
mass modification is by the substitution: 



T — > T c ff — VT 2 + (/V2tt) 2 . (1) 

where fib is baryochemical potential. 

A brief description of the thermal model particle spec- 
tra is given in the following section [2] An analysis of the 
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pion spectra based on the proposed model with medium 
mass modification is presented in section [3l followed by a 
summary and discussion of our work in section [4] 



flow [9|I10|. Thus the parameter T comprises, in a quali- 
tative manner, the effect of radial matter flow. Using in 
first approximation the Doppler formula, 



2 Particle spectra in the thermal model 

An important question governing the validity of the study 
presented here is: In what way is the pion spectrum mod- 
ified by the medium? When a medium-modified hadron 
emerges into medium free space, the hadron mass relaxes 
and the hadron picks up (or if appropriate releases) the en- 
ergy from/to the vacuum, and its mass returns to the nor- 
mal free space value. In such a process, within an isotropic 
medium, the translational momentum of a particle cannot 
change. However, along with the mass, the energy of the 
particle will also change. Since the momentum distribution 
of particles with p = E/v is not changed in the process of 
free space mass restoration: 

1) the number of thermal particles with momentum p is 
governed by the thermal momentum distribution of the 
source, obtained with medium modified mass, for exam- 
ple in the Boltzmann (classical) limit: 



d 3 N 
d 3 p 



-y/m(Ty+p*/T 



(2) 



and thus for particle spectra all masses are to be under- 
stood to be medium-modified; 

2) one should be able to reconstruct a decaying resonance 
using the medium modified masses for both the resonance 
and its decay products, for example: 

m\ (T) = [En + £ N ) 2 - (p* + Pn) 2 



ml(T) 



\(T) 



Ptt ■ Pn 



AT)- 



(3) 



Here all momenta are as measured in the final state. 

This discussion shows that in order to obtain the ob- 
served momentum spectra of produced particles, we have 
to use medium-modified masses in the usual thermal model 
expressions, and discuss the results as function of the mea- 
sured momentum. Other than the medium effect, our anal- 
ysis proceeds along previously established methods [9lfT0l 
lllj . The thermal momentum distribution of primary pi- 
ons is, by assumption, given by the standard thermal Bose 
distribution (here h = c = Ub = 1), see also Eq. ([2J: 



1 



dN 



2tt ptdptdy 



N 



Turrit cosh y 

77 g/3m t cosh y y 



(4) 



(3 = 1/T and N n is the normalization factor, containing 
among other effects, the volume of the source. Rapidity 
tanh - (p\/ E) and transverse mass mt = yrn 



2 +Pt 



are convenient covariant variables. In the center of the hot, 
dense collision zone, high pressure is produced. This pres- 
sure will cause a collective motion of the system, matter 



f 



f 



(5) 



for v ~ 0.25c and T t h — 100 MeV (a value generally ex- 
pected for the kinetic freeze-out in the reaction energy 
domain here considered) we see that the slope parameter 
of the spectra would be T = 129 MeV. For further details 
of emerging reference particle spectra we refer to Ref.[I2 . 

It is common to distinguish the pion charge polarity 
in the pion fugacity T v using the third component of the 
isospin ^3, (or alternatively the pion charge). The conven- 
tional definition of TV is [T3] : 



T 

■*- IT 



7 2 A/ 3 , 



T 



7 2 A 



g n I3 



T~± < e 



{3m 



(6) 



The last expression reminds us that is bounded by 
the Bose singularity. The presence of the parameter ^ q 
is required since it can influence the magnitude of pion 
charge asymmetry |llj when it approaches in magnitude 
the Bose singularity value. We recall that \i are particle 
fugacities and = e Mb / T . 

From Eq. (gj) we derive by integration over ptdpt the 
rapidity distribution: 



dN 
dy 



Turrit coshy 

oftmt cosh y y 



-Ptdpt- 



(7) 



We expect matter flow of the system in the longitudinal 
direction: 



dN(y) 
dy 



dr) 



dNjy - rj) 
dy 



f3 L = tanh(7y max ), (8) 



where f/ max = — ?y m i n , from symmetry about the center of 
mass, and (3l is the maximum longitudinal velocity. Thus 
we have the distribution of thermal pions 



1 dN(y) = ^ /"'" 
2ir ptdptdy "j^. 



Turrit cosh(y - rj) 



g/3m t cosh(y — rj) ^ 

Then the rapidity distribution is 

T^TOt cosh(y - 77) 



dr,. (9) 



dy Jr] 



o/3m t cosh(y-?7) _ y 



dptdr,. (10) 



3 Analysis of the data 

In order to model of the charged pion spectra [T] : 
1) We describe the high p t component of the pion mt spec- 
tra in terms of a thermal Bose shape, which in fact in this 
case is essentially an exponential. In this way we obtain 
the effective inverse slope T. This T comprises a combi- 
nation of the intrinsic thermal temperature T t h and the 
shift due to the collective longitudinal and radial expan- 
sion flow. The longitudinal flow is required for an accurate 
description of the rapidity spectra. 
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2) We form the difference between the direct pion spec- 
trum determined by the high m t fit and the experimental 
spectrum. We then model this difference in terms of the 
pion decay spectrum Z\ 33 (1232) — > N + n. We form the 
pion spectrum from direct and decay components using a 
fitted strength of the not-scattered Z^-decay component. 

3) We find that the pion spectrum can be described pre- 
cisely when we reduce all hadron masses using a uniform 
multiplicative factor. A further improvement is apparently 
possible when the radial flow is introduced, but we do 
not pursue such refinements here for reasons already dis- 
cussed. 

The thermal momentum distribution of primary pions 
is, by assumption, given by the thermal Bose distribution 
Eq. (H|). In the center of the hot, dense collision zone, high 
pressure is produced. This pressure will cause a collec- 
tive motion of the system which has been experimentally 
observed [T4JJT5] . However, it is not possible for us to dis- 
entangle thermal and collective flow effects p|]. Thus, the 
m^-slope we obtain comprises some Doppler-like shift due 
to source motion. Aside of the collective radial motion we 
expect that, in the longitudinal direction, the memory of 
the projectile-target motion will remain. 

Such longitudinal flow (3l is considered in Ref. pQ, em- 
ploying a flat longitudinal rapidity profile bounded by 
Vmax = -Vmin, where j3 L = ttmYi(ri max ). A somewhat sim- 
pler approach is to consider two fluids moving apart with 
a remainder of the original projectile and target rapidity: 



d 2 N(p t ,y) 
Ptdptdy 



l d 2 N( Pt ,y- 
2 ptdptdy 



l d 2 N{p t ,y + T] ) 
2 ptdptdy 



(11) 



We will fit 770 to the experimental rapidity distributions. 
For the mt spectrum inclusion of these two sources means 
that at y — we use Eq. ((4} substituting m t — > m t cosh^o- 

The parameters of the thermal model are presented 
in the top section of table [TJ where T s i opc = T/coshrjQ. 
T comprises the effect of transverse expansion dynamics. 
For large pt the pion fugacity is irrelevant: the thermal 
spectrum has a Boltzmann shape and is dominated by 
direct thermal component, and the fugacity becomes an- 
other component in the yield normalization. In this limit 
value of the mass of the pion is also irrelevant, as it is 
negligible compared to the value of p± . Hence the domain 
of large p t >> fixes the value of T s i ope . There re- 
mains a major difference between experimental spectrum 
and a single component thermal spectrum as was noted 
in pp. This difference is shown in Fig.[TJ This difference 
appears similar at all reaction energies and is localized in 
to a momentum range p t < 0.4 GeV. This suggests as the 
common and unaccounted mechanism of pion production 
the Zi33-resonance decay. 

For small p t there is considerable impact of the Bose 
nature of the pion on the spectra. Furthermore, the dif- 
ference in the shape between the positive and negative 
pion spectra can be in part accommodated by differences 
in Y„± , and in part by the associated difference in the 
relative yields of the zA 33 -resonances as indicated by the 
power of A/3 for each of the polarities, see last column in 
tabled The branching into the different channels for each 



Table 1. Top: statistical parameters used to describe 7v spec- 
tra; bottom: modified hadron masses as used at different beam 
energies. 



Beam Energy 


2A GeV 


4A GeV 


6A GeV 


SA GeV 


T slope [MeVj 


121 


128 


133 


135 


Vo 


0.15 


0.25 


0.3 


0.3 


T [MeV] 


122 


132 


139 


141 


A I3 


1.2 


1.15 


1.15 


1.15 




1.85 


1.62 


1.44 


1.44 


i^-[GeV" 3 ] 


120 


339 


521 


647 


i^ + [GeV- s ] 


136 


383 


587 


727 


TV^efB/Tch 


2,170 


983 


630 


536 


m/mo 


0.72 


0.63 


0.59 


0.57 


m A [GeV] 


0.880 


0.774 


0.722 


0.704 


m N [GeV] 


0.671 


0.590 


0.550 


0.537 


m^GeV] 


0.100 


0.088 


0.082 


0.080 


X7 dof 


1.16 


1.20 


1.04 


1.19 



of the ZA33 -resonances is relevant considering the isospin 
weights. Since we allowed for the isospin asymmetry in 
the fugacity, we could insist in normalizing the yields 
by same common factor. On the other hand, a consider- 
able improvement of the description of the spectra arises 
if the normalizations differ slightly, within the 15% uncer- 
tainty of the experimental data. The normalization factors 
we use are also presented in table [U Aside from experi- 
mental effects, such minor variation could also arise from 
differences in rescattering cross sections of charged pions 
in neutron-rich baryonic matter. 



Table 2. Delta decays and fugacities. 



Decay Channel 


r A 


A++ -> (tt++p) 
A+ -» l/3(7r++n), 2/3(tt° + p) 
A -» 1/3(tt- +p), 2/3(7r° + n) 
A~ — (77" + n) 


ll^ 
3,-1/2 

!q A /3 



To describe this difference spectrum quantitatively, we 
consider ZA33 two body decay, and employ the well known 
results, see, e.g., Ref. [91H0] . 



dN, 



_ m R b 

Ptdptdy inp* 



dyi 



mjCOsh 2 (y 



2( + ) 



dm 2 tR 



vr) ~ 

dN R 



pj 



, 2 (-> 

tR 



( m tR ~ m tR)(m tR 



f-K m t Rdm t Rdy R 
"Hr ) 



,(12) 



where b is the branching ratio, and p* is the momentum 
of the decay particle in the Z\33-rest-frame. For the ZI33- 
resonance spectrum we assume at first a Boltzmann dis- 
tribution, with T determined by the pion spectrum and 
the Zi33-fugacity shown in table [2] 



1 



dN, 



2tt m tR dm tR dy 



= N A e 



T A m t Rcosh y R 
exp(/3m t _RCOsh y R ) 



(13) 
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-2 1 1 1 1 1 1 1 -2 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

p, (GeV) p t (GeV) 

Fig. 1. Difference of experimental tt^ yields and the direct thermal (th) yieids, for the mid-rapidity data bin pQ. The contribution 
from mass modified A33 decay is shown by solid line. Contribution from mass modified A33 decay with radial flow is shown by 
dashed line for the case of 8 GeV beam energy (top panels). Dot-dashed line shows the A33 decay contribution without mass 
modification. 



In the resulting pion decay component we also introduce 
the longitudinal flow according to Eq. (JTTJ). In this way we 
allow effectively for the longitudinal motion of the ^33- 
source. The total pion distribution is composed of the di- 
rect and the decay contributions. 

When computed with the free space masses of tt, ^33 , N 
the resulting spectrum of decay pions is found to have 
the shape seen in Fig.[IJ but is systematically shifted to 
a higher momentum than the experimental data. In order 
to reduce the energy scale to the level seen in the differ- 
ence spectrum in Fig.[TJ we scale down the three hadron 
masses 7r, ^33,^ by the same factor, which is chosen in 
a qualitative manner. The resulting masses we use are 
stated in the bottom section of table [TJ the modification 
is non-negligible. There is a tendency for the mass reduc- 
tion effect to increase with collision energy. As is seen in 
Fig. [H (solid lines), after this modification we describe the 
experimental data rather well. 



We also have to allow for the presence of matter flow. 
In general the radial flow has the effect of flattening the 
^33 spectrum more than the tt spectrum, since the ef- 
fective inverse slope T, which combines the intrinsic tem- 
perature T t h with the blue-shift of flow, grows somewhat 
with particle mass and collective velocity v 16J. We ac- 
count for the difference that the flow effect has on tt, and 
the much heavier ZI33 as follows: We keep all quantities 
except for the distribution Na unchanged. We modify the 
^33-spectrum using a simplification of the model proposed 
in Ref. [17]. For the flow rapidity profile rj t (r) — rj t f(r/ R) a , 
we take a = and -qtf = 0.3. We allow here a cylindrical 
fireball since this is not a significant element in our ap- 
proach and the model is simple and transparent. The con- 
tributions from the Z\33 decay with radial flow are shown 
in Fig.[T] solely for 8 GeV beam energy, where this effect 
can be expected to be largest, as dashed lines. There is 
a mild improvement in the high p t part of the difference 
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spectrum. These results suggest that the complete and 
detailed inclusion of the radial flow, which has the largest 
influence for highest reaction energy, will further improve 
the understanding of the difference spectra seen in Fig.[TJ 



The integral of the transverse momentum spectra leads 
to the rapidity distribution shown in Fig.[2j In perform- 
ing this integral it is relevant for the rapidity dependent 
normalization to disentangle in the slope temperature the 
flow effect from the true temperature. Since we already 
have seen that matter flow is desired, it is not surprising 
that when we do this introducing for simplicity a longi- 
tudinal flow tjq, the rapidity spectrum is described very 
well. The resulting % 2 /dof is given in the bottom of table 
[TJ It should be remembered that this is not a best fit op- 
timizing the 8 parameters (4 statistical, 3 normalization 
parameters, and one mass scale parameter which deter- 
mines the three hadron masses) at each reaction energy, 
and that we have treated the flow in somewhat cavalier, 
but as we think, precise enough manner. 



4 Discussion and Summary 

We have obtained an accurate description of pion 7r ± spec- 
tra, consistent across the four reaction energies considered 
here. Our objective we had here was to show that the pion 
spectra we consider are sensitive to the introduction of the 
mass change parameter. We believe that we have demon- 
strated this in that we can explain quite satisfactorily the 
spectra considered. The relevant model parameters were 
determined in part by using the common knowledge and 
in part optimizing the parameters to the data. We did not 
make an effort to minimize the value of x 2 j since several 
important physics elements are beyond the scope of the 
current approach. 

We generally proceeded as follows: we chose a reason- 
able value of the hadron mass reduction parameter, the 
longitudinal flow 770, and the two fugacities. Once this 
choice is made, the temperature T (common for both tt + 
and 7r~ ) , the two pion normalization factors N^± and the 
^33-yield normalization are fitted to the data. We then 
tried another reasonable set of the three initial param- 
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eters in an effort to further reduce the \ 2 an d did this 
several times. 

We would like to make a few additional remarks about 
the physical meaning of the parameters, possible future 
directions, and consistency of our approach: 

1) One could argue that the softening of the pion decay 
spectrum, i.e., the difference spectrum seen in Fig. Q] is not 
due to mass modification but is due to decay pion scat- 
tering in matter. It is however hard to understand why 
this process at all reaction energies would have the effect 
of producing the Z\33-decay spectrum shape, with reduced 
energy scale. We note that though the data description we 
present is relatively good, at low pt there is some poten- 
tial to improve the model, e.g. by considering the effect of 
pion scattering. 

2) In table [TJ Na is common for different isospin states 
of A For 6, 8 GeV beam energy cases, Na and N v are 
comparable. However, Na > N n for 2 GeV beam energy. 
This suggests that the direct thermal pion component is 
dominated by the indirect pion production via the A33- 
resonance. Indeed, as is seen in Eq. (fT"3f . the ^33-yield 
comprises a factor, e tlB ^ Tch , where /xb is the baryon chem- 
ical potential. The rapid growth in the normalization fac- 
tor with decreasing reaction energy is due to the expected 
increase of /is /Teh- 

3) Hadron model calculations [5] find mass shifts which 
are small, and decay width increases which are large. We 
note that the effect obtained in the bag model effect dif- 
fers. The decrease in the hadron masses is accompanied 
by a decrease in transition energy, and thus the decay 
phase space diminishes. This reduces the width. There is 
a further potential reduction due to the reduction of the 
intrinsic decay matrix element, since for a smaller value of 
the bag constant B the quark bag grows in size, diluting 
the strength of the reaction matrix element. We recall that 
at present an increase of mass shift has not been observed 
experimentally while a decrease in width is not inconsis- 
tent with the available data. Thus a priori, a bag model 
description seems to agree better with data. 

4) The consistency of our model can be tested by intro- 
ducing the actual /^-resonance spectra obtained from a 
reconstruction by the invariant mass method of the mass 
modified /^-resonance using modified masses of decay 
pions and nucleons. Success of this would indeed be a 
convincing step supporting the path of analysis proposed 
here. Also, it would allow us to include in our approach 
the spectra of nucleons. 

5) The matter flow field arises from initial conditions and 
ensuing dynamics driven by equations of state. Thus there 
is in principle a relation between the magnitude of the flow 
and the freeze-out condition (temperature, matter den- 
sity) for each heavy ion reaction energy, if initial condi- 
tions could be understood. Our results are consistent with 
what may be called, in absence of theoretical understand- 
ing, the general wisdom about the magnitude of matter 
flow, temperature and density expected in conditions con- 
sidered here. 



In summary, we have studied the precision spectra of 
7r obtained for kinetic beam energy per nucleon 2, 4, 6, 8 
A GeV at AGS. We can interpret the it spectra as origi- 
nating in a thermal fireball of dense matter. The necessary 
ingredients are aside of direct tt^ 1 production, secondary 
Z\33-resonance decay, longitudinal flow, and, a significant 
mass modification of the involved hadrons. In fact, the 
mass modification is the new physics element in the anal- 
ysis presented, allowing a significant improvement of the 
data description using a single thermal source. 
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